Global vs local Casimir effect 
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Abstract 



This paper continues the investigation of the Casimir effect with the 
use of the algebraic formulation of quantum field theory in the initial value 
setting. 

Basing on earlier papers by one of us (AH) we approximate the Dirich- 
let and Neumann boundary conditions by simple interaction models whose 
nonlocality in physical space is under strict control, but which at the same 
time are admissible from the point of view of algebraic restrictions imposed 
on models in the context of Casimir backreaction. The geometrical setting 
is that of the original parallel plates. By scaling our models and taking ap- 
propriate limit we approach the sharp boundary conditions in the limit. The 
global force is analyzed in that limit. One finds in Neumann case that al- 
though the sharp boundary interaction is recovered in the norm resolvent 
sense for each model considered, the total force per area depends substan- 
^j , tially on its choice and diverges in the sharp boundary conditions limit. On 

the other hand the local energy density outside the interaction region, which 
in the limit includes any compact set outside the strict position of the plates, 
has a universal limit corresponding to sharp conditions. This is what one 
should expect in general, and the lack of this discrepancy in Dirichlet case is 
rather accidental. Our discussion pins down its precise origin: the difference 
in the order in which scaling limit and integration over the whole space is 
carried out. 
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1 Introduction and the main idea 

The most natural setting for the consideration of the Casimir effect is the algebraic 
approach. This approach allows a mathematically rigorous analysis of the effect 
and gives a clear understanding of the sources of the difficulties one encounters in 
more traditional treatments. In application to quantum fields this analysis rests, 
in broad terms, on the following cornerstones. 

(i) A quantum relativistic theory is defined by an algebra of observables, in simple 
cases defined directly by 'fields' (scalar, electromagnetic). 

(ii) Each particular physical system obeying this theory is described by a Hilbert 
space representation of this algebra. Inequivalcnt representations refer to physi- 
cally non-comparable systems or idealizations (such as a local isolated system and 
a thermodynamic limit system). 

(iii) The change of external conditions under which a quantum system is placed 
leads to a change of the state of the system. 

For the calculation of the global Casimir-type effects, as the backreaction 
force, one needs models which respect the above three constituents of a quantum 
theory. Thus the model of a quantum field should be based on one definite algebra, 
and the interaction with the external conditions should not lead to a change of 
its representation. If this condition is fulfilled then the Casimir force results from 
the change in the expectation value of one and the same energy observable, as 
defined by free field, as the state changes with changing external conditions (such 
as position of macroscopic bodies) . 

This analysis has been conducted at length by one of us in [I] and [2] , where 
also clear cut criterions for the admissibility of external interaction models for 
a class of systems were formulated. In application to the free quantum scalar field 
these amount to the following. Let cf>(x) be a scalar field and denote by h 2 the 
standard selfadjoint extension in L 2 (]R 3 ,d 3 x) of —A. The free field dynamics is 
then 

(d 2 + h 2 )^(t, x) = o . (i) 

Suppose now that one introduces to the system external macroscopic bodies, such 
as conducting plates in the original Casimir system, which change the dynamics 
of the field. Let a denote free parameters of these bodies, such as separation of the 
plates, and let the modified dynamics (for fixed a) be given by 

(d 2 + h 2 a W(t,x)=0, (2) 

where h 2 is a positive selfadjoint operator. Now, the two settings can be described 
by one choice of observables algebra, and in common representation of this algebra 
if, and only if, h a (h a — K)hr x l 2 is a Hilbert-Schmidt operator in L 2 (R 3 , d 3 x) , 
that is 



Tr 



h^^iha-h^iha-h)^ 1 ^ 2 <0O. (3) 



Suppose this condition holds and let the algebra of the field be represented in 
some Hilbert space T-L. Let further H be the energy operator as defined by free 



field dynamics, and let fl a be the minimal energy state vector as defined by the 
modified dynamics ([2]) . The Casimir energy is then given by 



£ a =(Q a ,Hn a ) = jTr 



(ha-tyhzHha-h) 



(4) 



That this energy be finite is another condition on the model of h a , and only if both 
conditions are satisfied the Casimir problem has a finite solution and the Casimir 
force is then 

da 

These admissibility conditions say, roughly, that the modified dynamics h a 
cannot differ much from the free dynamics h 2 . Introducing sharp boundary con- 
ditions, such as Dirichlet/Neumann conditions on plates, violates these demands. 
One faces therefore the problem of an appropriate approximation for the descrip- 
tion of such plates. 

We consider the simplest geometrical situation, the original Casimir problem 
of two infinite, parallel plane plates at a distance a from each other. We assume 
that z-axis is perpendicular to the planes and the modification of dynamics affects 
this direction only, thus 



h* = ht + hl, h\ = h\ 



h 2 



where h\ is the free dynamics in the directions perpendicular to z-axis and h 2 a 
is a modification of h 2 = —d 2 in L 2 (R, dz); we refer the reader for details to [2]. 
In this setting the conditions of finiteness of (|3|) and (|4|) cannot be expected to 
hold as they stand because of translation symmetry in the planes, and must be 
replaced by conditions 'per unit area' of the planes. It has been shown in [2] that 
this amounts to 



Tr 



(h za - h z f 



< oo . 



Tr (h za - h z )h z (h za - h z ) < oo 



(5) 
(6) 



where the trace refers to the Hubert space L 2 (1R, dz). If these conditions are sat- 
isfied, the energy per unit area is finite and reads 



e n = Tr 

a 24tt 



(h za ~ h z )(2h z + h za )(h za - h z ) 



(7) 



A class of models for h 2 a imitating the boundary conditions, but con- 
sistent with the above demands was considered in [2\ The idea was to take 
h 2 a — (function of)(h z ,hf a ), where (/if a ) 2 is —d 2 with boundary conditions at 
z = ±a/2. The choice of functions assured that for small spectral values of h 2 
and {hf a ) 2 the models reproduced the sharp boundaries, while for large spectral 
values tended to free dynamics. Moreover, one could introduce a scaling parameter 



fi such that for \x — > oo the models approached the sharp boundaries in the whole 
spectrum. For the Casimir energy per area in the rescaled models e^ (we prefer to 
work with A = 1/fJ, here) one then found 

£ ^^ + ^- Ii ^ + (tcrms^0forA^0), (8) 

where £00 and c are constants, c = in Dirichlet case, but c 7^ in Neumann case. 
It was also shown that the direct sum of the two models describes the setting of 
the electromagnetic field between conducting plates. Thus the Casimir force per 
unit area is then 



d[e*(D) + e*(N)} 2c 



da Aa 3 240a 4 



(terms — >• for A — >• 0) . 



The second term reproduces the well-known Casimir's formula, but the first term 
is model-dependent and dominates for large a. Moreover, in typical situations there 
is c > and the force becomes repulsive for large a. 

In the present paper we want to find out whether these results will be con- 
firmed in another class of models, constructed in a wholly different way. Rather 
than manipulate spectral properties directly, now we want to approximate inter- 
action with the plates directly in the physical space. It is easy to see that strictly 
local potential interaction of the form [h1 a ip](z) = —d^ip(z)+V(z)tp(z) violates our 
conditions. Therefore we replace V by a slightly nonlocal integral quasi-potential 
[V^](z) = / V(z, z')ip(z')dz' with the kernel V(z, z') concentrated around the po- 
sition of the plates. We show that with an appropriately defined scaling of some 
simple kernels of this kind one can reproduce sharp boundary conditions on the 
plates in the limit. The Casimir energy can again be calculated and for a class of 
models the result (jHJ is confirmed. However, in general the Neumann case proves 
to be even more singular here than in the models considered in [2] and the uni- 
versal term could be disturbed. Nonlocal quasi-potentials has been considered in 
the Casimir context by other authors before, but in different formalisms and with 
rather different motivations (see e.g. [3]). 

The present choice of models makes also possible a local analysis of the local 
energy density. We show that outside the interaction region the density tends in 
the scaling limit to a well defined universal form corresponding to sharp boundary 
conditions. The present mathematically rigorous setting allows the comparison 
and better understanding of the local - global relation. The model-dependent and 
divergent (in the limit) contributions to the global force are due to the interaction 
region. We discuss this point more fully in the Discussion section. 

For a more extensive discussion of the background of the present paper, as 
well as for more extensive literature we refer the reader to pQ and [2J. We define 
our models in Section [2] Appropriate scaling of these models is shown to reproduce 
the sharp boundary conditions in Section [3] Spectral properties of the models are 
discussed in Section |4] and the admissibility of the models in the sense mentioned 



above is proved in Section [5j It is shown in Section [5] that the Casimir energy of 
the scaled models is obtained by the expansion of the formula for energy in inverse 
powers of a, and this expansion (up to a significant order) is obtained in Section[?l 
For comparison, in Section [5] we obtain local results and their scaling limit. The 
discussion occupies Section |H1 More technical points of our derivations are shifted 
to Appendices. 

2 The models 

We postulate for our analysis the following quasi-potentials 

V = a(\U b g){U b g\ + \U- b g){U- b g\), b = a/2 > , a = ±1 , (9) 

where U is the translation operator and a = 1,-1 corresponds to Dirichlet (D) 
and Neumann (N) conditions respectively. These conditions will be achieved in 
the two cases by an appropriate scaling limit to be defined below. In all what 
follows one should keep in mind that unless stated otherwise we treat parallelly 
both cases, but the dependence of quantities on a is suppressed. 

In position representation the quasi-potential is an integral operator 
(Vip)(z) = J V(z, z')ip(z') dz' with the kernel 



V(z, z') = g [g(z - b)g(z> - b) + g(z + b)g(z' + b) 
For the functions g we assume that 

[-*&/(*) lfcr = - 1 : (N) 

where / is a complex, compactly supported smooth function, with the following 
properties 

f(-z) = f(z), 8uppfC(-R,R),R<b, /(0)^0, (11) 

11/11=1, (N) (12) 

where / is the Fourier-transformed function 

/(p) = -L ff(z)e-^'dz, 

V 2vr J 

and the last property is assumed only in the Neumann case. The first condition 
reflects the symmetry of each of the plates, the second says that the nonlocalitics 
of the two interaction centers at z — ±6 do not overlap, and the third and fourth 
arc technical. 



We denote by h z and h za the selfadjoint, nonnegative square roots of the 
operators 

h 2 z =-(d/dz) 2 , h 2 za = h 2 z + V, (13) 

respectively. Operator h 2 is the standard one-dimensional Laplace operator (with 
opposite sign) , while h 2 a is its Kato-Rellich perturbation, with unchanged domain, 
as V is bounded. The (strict) positivity of h 2 a in the Dirichlet case is obvious, while 
in the Neumann case one has by (fTUf that for all ip in the domain of h 2 there is 

OA, h 2 za i>) = ll^'H 2 - \(u +b f^')\ 2 - |(c/_ b /, V')l 2 > (N) 

where ip'(z) — dtp{z)/dz. The functions U± b f do not overlap, and satisfy 
ll^ifc/ll = 11/11- Thus by (fT2|) they form an orthonormal system, which implies 
(ip, h 2 a ip) > for each ip. (Here a more strongly bounding condition ||/|| < 1 
would produce a strictly positive operator; this, however, would not lead to the 
recovery of Neumann condition in a limit to be defined below.) 
For w 2 £ <D with 3m w 2 ^ the resolvents denoted by 

G (w 2 ) = (w 2 - hi)' 1 , G{w 2 ) = (w 2 h 2 za )- 1 , (14) 

are bounded operators. In all what follows for given w 2 we fix w by 3mw > 0. We 
also introduce the T-operator known from the stationary scattering theory 

T(w 2 ) = V + VG Q (w 2 )T(w 2 ). (15) 

This equation may be explicitly solved for T: making the Ansatz 

TV) = ( \U +b g) \U_ b g) ) T(w 2 ) ( J^ 9 ' ] , (16) 

\{U-b9\) 

with T(w 2 ) a numerical matrix, one easily finds 

fa-( gi Go(w 2 )g) -(U a g, G (w 2 )g)\ _1 
T(w 2 ) = . (17) 

V-(C/ Q .g,G ( W 2 ) 5 ) a-(g,G (w 2 )g)J 

The resolvent may be then expressed by 

G(w 2 ) = G (u; 2 ) + Go(u; 2 ) 7> 2 )G (™ 2 ) . (18) 

In momentum representation, taking into account (jlQI) and denoting 

?v = P^/iW (e" lbp e+^) , M p = \f(p)\ 2 , (19) 



we have 



(p\T(w 2 )\q)=T p T(w 2 )?l, (20) 

<p\ G(w 2 ) - G„( W 2 ) | ? ) = —5L^T( W 2 )^^ , (21) 

w — p w — q 



with elements in the matrix (|17p given by 

/iap \-g 71, r 
^-f ^dp = nrar<V a *'M„ ) . (23) 

The integral in the last formula is calculated in the complex plane by residues, with 
the use of analyticity and asymptotic properties of M w discussed at the beginning 
of Appendix |XJ 

3 Reproduction of the sharp boundary conditions 

We consider now a family of rescaled quasi-potentials V\, A 6 (0, 1), built as in 
(|9|), but with the use of rescaled functions g \ instead of g. We write the scaling in 
several equivalent forms: 



fix 



(z) = X'ig (|) , f x (z) = A" 1 "*/ (|) , f x {p) = A-*/(Ap) , 



and note also that M Pt \ = \~ a M p . The rescaled potentials give rise to the cor- 
responding operators h za ,\- All quantities referring to these operators acquire the 
subscript A. 

Let h B a be the selfadjoint, positive square root of the operator 



d 2 

,S 12 _ " 

dz 2 



\h B \ 2 = 

V l za\ 



(24) 

boundary conditions 



with standard domains in L 2 (R), Dirichlet (for a = +1)/Neumann (for a = — 1) 
conditions in z = ±b, and denote by G B (w 2 ) the resolvent of [h B a ] 2 . Our objective 
in this section is to show the limiting property 

lim \\G x (w 2 )-G B (w 2 )Us = 0, (25) 

where the Hilbert-Schmidt norm is ||A||h S = Tr(A*A) > \\A\\ 2 . From this relation 
it follows then the norm convergence 

lim jF(h za , x )- F(h* a ) ||=0 

for each continuous and vanishing in infinity complex function F on 1R, and the 
strong convergence 

]im\\[F(h za>x )-F(h* a M\=0, 

A^0+ 

for each bounded continuous function F and vector ip G L 2 . 



It is clear from the form of Eq. (|2ip that G(w 2 ) — Go(w 2 ) and its scaled 
version G\(w 2 ) — Gq(w 2 ) are finite rank, hence Hilbert-Schmidt, operators. Thus 
it is sufficient to calculate the strong-L 2 (]R 2 , dp dq) limit for A — > + of the integral 
kernel (p\ G\(w 2 ) — Go(w 2 ) \q). First we consider the numerical matrix 7a, and for 
later use we also look at higher orders in A. We observe that 



M 



x P 



Mn 



dp = I + 0(A) , 



with 



h = 



M - M, 



P 



^dp. 



(26) 

(27) 



This is shown by writing the difference of (|26j) and (|2T[) as 

M u - M Xw 



Xw 



(Xw) 



■ du 



and Fourier-transforming the integral as a scalar product of two L 2 -functions. This 
integral is shown in this way to be bounded by a constant. Using (|2"2"1) . ([2"5|) . ([2"B)l 
and the assumption (JTSJ) (in Neumann case) we find 

a - ( 5A , G (w 2 )g x ) = tTr(Xw)-°M + (Xw) 1 ^ (^ - I ) + 0(X 2 -°) , (28) 

- (U a g x , G (w 2 )g x ) = iv(Xw)-° e iaw M + 0(X 2 -°) . (29) 

From these we get 



Tx(w 2 ) = 



-i(Xw) a 



ttM (1 



a 2iaw 



+ (Aw) 1+(T (i±2. _ J ) x {matrix independent of A} + 0{X 2+a ) 



Next, we observe that (this is shown in Appendix [AT)) 



^° 2 h(p)p~ 



f($)p- 



w — p w — p 

and in addition, for Neumann case 

A-*/a(p) 7(0) 



< 



{const (w) A 2 
const(ui) A3 



(D) 

(N) 



w 2 — p 2 



w 2 — p 2 



<const(u;)A5 (N) 



(30) 



(31) 



(norms of functions of p as elements of L 2 ). Now we easily obtain in the 
L 2 (K 2 , dp dq)-sense 



lim (p\Gx(w 2 )-G Q (w 2 )\q) = 



iw 



(pqY 



A^0+ 



7r(l - e 2iaw ) {w 2 - p 2 )(w 2 - q 2 ) 



e -ibq e ibp _ e iaw e ibq e ibp _^_ p ibq & -ibp _ ^iaw e ~ibq & -ibp 



We transform this to position representation and get 



s- lim (z\G x (w 2 )-G (w 2 )\z') 



2iw (1 - e 2law ) 
x 



6(b + z)e i{ - b+z ^ w + a9{-b - z)e~ i( - b+z '> w 

aO{-b - z')e- i{b+z ' )w + 9(b + z')e l{b+z '^ w 
-<r6(b- z') e iaw e ^ b ~ z ') w - e(-b + z 'y^ e -i(b-z')w 

6{-b + z)e- l{b - z ^ w + ad(b - z)e l{b - z ^ w 

< \aO{b - z')e i(b - z ' )w + 9(-b + z') e - i(b - z ' )w 

- a 6(-b - 2 ') e «™ e -*(<>+2> _ Q(b + z 'y»w e i(b+z')w 



We shall use also the explicit form of the unperturbed Green function in this 
representation 



(z\ G {w 2 ) \z') = - — 9{z ~ z')e l{z - z ' )w + 6(z' - z)e i{z '- z)x 
In this way we find 



(32) 



s - lim (z\ G x (w 2 ) \z') = (z\ G B {w 2 ) \z') , 

A->0+ 



where 

(z\G B (w 2 )\z<) 



(z\G Q (w 2 )\z'} + *-(*+*'+<>)- 
2w 



X(-oo-b){z)X(-oo,-b){z') 



(z\G (w 2 )\z')+a^ z+z '-^ w 
2w 



X(b,+oc)(z)X(b,+oo)(z') 



(z\G Q (w 2 )\z') + 



, ., ,. i { cos(zw) cos(z'w) sin(zu>) sinfVtu) 



w V 1 + ae 



1 — ae 



X(-b,b){z)x(-bM)(z') 



(33) 



and xn 1S the characteristic function of the set O. 

In the three regions (z\ G B (w 2 ) \z') differs from (z\ G (w 2 ) \z') only by solu- 
tions of homogeneous equation and satisfies the boundary conditions 



(±b\G B (w 2 )\z') = (z\G B (w 2 )\±b)=0, 
4 (z\ G B (w 2 ) \z') \ z=±b = £ (z\ G B (w 2 ) \z') |,,=±6 - , 



dz 



dz' 



(D) 
(N) 



so it is indeed the Green function of the Dirichlet /Neumann operator and therefore 
(|2"5|) is finally proven. 



We now want to acquire some information on the rate at which the limit ([55 
is achieved. Using (f30|) one finds for any p, r\ <G L 2 



(34) 



( r , M / rB( 2 , . J(l-/o)0(A) + 0(A 3 / 2 ), (D) 
{if, G x (w )rj) = {<p, G (w )r,) + i /2) ^ ^ 



The Neumann case turns out to be here, as in many other problems, more singular. 
However, we also note that if we assume that tp and r\ are in the domain of h z 
then the estimate (|31l) implies 

(<p, G x (w 2 ) v ) = (ip, G B (w 2 ) v ) + I O(X) + 0(A 3 / 2 ) . (N) (35) 

In the following two sections we treat unsealed models. The scaling is again 
considered in Section [BJ 

4 Spectral analysis 

We add now some further assumptions on the choice of functions /. We denote for 

fceiR 

M *- M >, (36) 



p 2 — k 2 
and demand that 

< I k for k ^ (D,N) , (37) 

h < 1 (D) . (38) 

Note that by continuity 1$ > (this is the quantity introduced in (|27|l). We also 
denote 

«N h = w{ a + / ^-_p% } = la + *)i*r - 1*1/* • (39) 

The operators /i 2 a are nonnegative, and outside a compact set in R they 
act as — d 2 . Therefore their continuous spectrum covers the whole positive axis 
and thus the spectrum is (0, +oo). This does not resolve the question of point 
spectrum, and we treat it first. 

The eigenvector equation 

/i 2 >fc = fc 2 V-fc, Vfe€i 2 (R), fc>0 (40) 

is solved in momentum space. It is easily seen that the distributional solution 
which is square-integrable at infinity must have the form 

~ ri,p~ ib P 4- r uP ib P i „ ^ 

mp)= b /_-: p^m, (4i) 

10 



with constants c±b to be determined. Putting this form back into Eq. (14U1) one 
finds for fc > that the constants c±b have to satisfy the linear system 

N k c +b - M k s'm(ak) c_ b = , M k sin(afc) c +b - N k c_ b = , (42) 

where the integration leading to coefficients M k sin(afc) is carried out with the 
use of analyticity and asymptotic behavior of M k discussed in Appendix [S] Now, 
for fc > the conditions (j3"T)l and (1551) imply N k ^ 0, thus nontrivial solutions 
to the system (|4"2"|) exist only if M k sin(afc) = ±N k , and in that cases (|4"Tj) takes, 
respectively, the form 

^ e ib P zp e~ ibp i-o- — 

V>fcb) = C-6 p2 ^ fc2 P~I(P) ■ 

The condition ^ fc e L 2 (R) requires that f(±k) = or e i6fc ± e" lf,fc = 0. Each of 
these cases implies M k sin(afc) = 0, which shows that there are no eigenvectors for 
fc>0. 

For fc = in Dirichlet case the solution (|4Tj) cannot be in L 2 as /(0) ^ 0. For 
Neumann case one finds that (|41[) is a distributional solution for any constants 
c-ti,, so we are free to choose them so as to satisfy the square-integrability of (|4"Tj) . 
This happens only for c_ b = — q, and then 

$ (p)=Af S J^Mf(p), (N) (43) 

where J\f is a proportionality factor. The normalization condition 1 1 ijjq \ \ L 2 = 1 gives 

W? = J T ■ (44) 

Summarizing, there are no bound states for Dirichlet case, however for Neumann 
case there is one bound state, which corresponds to the zero eigenvalue, described 
by (03J and flU}. 

We now consider the continuous spectrum and for this purpose use the sta- 
tionary scattering formalism. The improper eigenfunctions of scattering states in 
momentum representation are given in standard notation by 

r n , ,,,* x , n , (p\T(k 2 + tO)\k) 

ipk{p) = (p | k+) = <$(p - fc) + fc2 _ 2+ij() , (45) 

where T(w 2 ) is the operator discussed in Section [2 The variable fc takes all values 
fc 7^ and each spectrum point fc 2 has two-fold degeneracy corresponding to ±fc. 
Taking into account the results of Section [5] we can write 

\ k \o (M k -iN k M k e ia \ k \\~ X 
Tik 2 + i0) = 4M . (46) 

i7r \ M t e m l fe l M h -iN h 



ff 



For later use we write this in two alternative forms. We introduce matrix notation 

(47) 



7T 



M( P ) 
p 2 — k 2 



^ + V I ~^h d P 



.,1.1-.,, / 1 cos(afc) 



m =M(k) + M(-k)=2k^M k ^ o ; {ak) -pj. (48) 

and then 

r(k 2 +iO) = ^(M(k)+iC{k)Y l , (49) 

|(p|T(fc 2 +iO)|fc)| 2 = Tr[X(p)T(fc 2 +iO)X(fc)T(fc 2 +iO) t ] . (50) 

Finally, we calculate the inverse in (|46j) and write the result in the form 



nf-HQ-j- . ' tr \,J ' ,, " <fte ""V (5i) 

ITT 1- (g fc e M l fe l) 2 y_ gfce »|fe| ]_ y 

where 

1 M k 

Sk = i Qk = — MhSh ■ 

M k - iN k ' qk M k - iN k 
Some properties of s k function are shown in Appendix El 

5 Hilbert- Schmidt properties 

In this section we show that our model satisfies the admissibility conditions |5J 
and ([6]). If we write TR T for the l.h.s. of these two conditions, with r = for ([5]) 
and r = 1 for ([5]), then 

TR T = J \p\ T \ (p\ h za - h z \k+) \ 2 dkdp+^ J \p\ 2+T \Mp)\ 2 d P , (52) 

R 2 R 

where the first term results from the continuous spectrum space of h za and the 
second is the bound state contribution in Neumann case. The second term is 
evidently finite (by (|43jl ). and we restrict attention to the first one, which we 
denote TR™ n . In momentum representation we have 

, ,, \ mi n? / x (p\T(k 2 +iO)\k) 

(p| h za - h z \k+) = (\k\ - \ P \)Mv) = \ p \ + \k\ ' 

thus by change of variables for negative arguments we get 

TR ™ nt = / T-T1A2 E I (±Pl T ( fc2 + i0 ) \ ±k ) I'* d P < ( 53 ) 

J (P + K) ±± 



12 



where the signs in 'bra' and 'kef are uncorrelated and the sum is over all four 
possibilities. From now on we assume that k,p > 0. Using Eqs. (j48|) . (|50|) we can 

write 

Y, | (±p| T(k 2 + *0) |±fc) | 2 - Tr [£(p)T(fc 2 + iO)£(k)T(k 2 + i0)^] , 
±± 

and substituting here (from (|39|l ) 



£(fc) = _ (T(fc 2 + iojt- 1 - T(fc 2 + »o)- 



we give Eq. ([53]) the form 



rp-pcont 



/cp T 



T J (p + ^) 2 

HI 



mc 



iTr [£(p)T(fc 2 + i0)] dfcdp 



(54) 



With the use of 



kme 



iTi [C(p)T(k 2 + i0)] 



and (|5ip we have 
4 



pi-^fci+^Mp SKe 



Sft 



1 — cos(ap)qke lak 



1 - (tffce 



iak \2 



(55) 



Writing out the real part gives us the appropriate behavior of the nominator for 
k = and the whole expression becomes proportional to M&. Using the estimates 
(|53|) . (HU and <JHZ1) one finds that TR™ nt are finite and the admissibility conditions 
([5]) and ([6]) are satisfied. 

6 Scaling 

We return to the scaling transformation to view it from a different point. If we 
make the dependence of the potential V on a explicit by writing it as V a and the 
rescaled potential as V a ,\ then we have 

V a , x (z,z')=\- 3 V a/x (z/\,z'/\). 

It is then an easy exercise to show that this implies a simple scaling law of the 
eigenfunctions (|4"3")) and (1431) : 

#£ 6 (p) = \$ Xk , b/x (\p) , $£ 6 (p) = X 1/2 ^ b /x(Xp) 

(different powers of A reflect different normalizations: to the Dirac delta and to 
unity respectively). Denoting the scaled versions of (|52|) . with explicit dependence 
on a, by TR^ a we find 

TR T ,a = X T TR Ta / A . 
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7T 2 


487rA/ a 2 87r 2 M 2 a 3 


1440a 3 



Thus the admissibility conditions are satisfied also for the rescaled potential (but 
not in the limit). In the same way one obtains the scaling law for the Casimir 
energy: 

£a=A~ 3 e aA . (56) 

Therefore to identify the scaling behavior of the energy in the limit it is sufficient 
to expand the unsealed energy in inverse powers of a up to the third order. This 
will be done in the next section. 

7 The energy 

In this section we prove the following expansion of Casimir energy for large a: 

(D) (57) 

+ o(a- 3 ) . (N) (58) 

We postpone the discussion of this result to the concluding section and here only 
note that there are functions in our class for which Iq = 0, and then the a~ 3 -term 
has the known universal form. 

The two conditions ((5]) and (j6]), considered in Sectional imply already finite- 
ness of the Casimir energy per unit area ([7]), as mentioned in the Introduction. This 
can be easily seen: we observe that conditions ([5]) and ([6]) mean that A = h za — h z 

1/2 

and h z A are Hilbert-Schmidt operators. Also, from 

Ah za A = Ah z A + A 3 

1/2 

we infer that h za A is HS as well, which is sufficient for the claim. 

The expression ([7]) closely parallels that of the condition ([6]) and can be 
written in analogy to (f5"12j). We split the trace in ([7]) into two terms and calculate 
these traces in h za or h z (improper) basis in the first and second line below, 
respectively. Then we insert the spectral decomposition of the operators h z or h za , 
respectively. In this way we get 

Tr [2Ah z A] =2 J \p\ \ (p| h za - h z \k+) \ 2 dkdp+(l - a) J \p\ 3 \Mp)\ 2 dp , 

R 2 R 

Tr [Ah za A] = f\k\\ (p\ h za - h z \k+) \ 2 dkdp . 

R 2 

Therefore the expression ((Tj) may be written in the form 
£a = i / (2bl + |fc|) ' {PI Hza " hz lk+) \ 2dkdp+ ^~ J \P\ 3 \Mp)\ 2 dp- (59) 

R 2 R 
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The bound state contribution e^° und - the second term - will be calculated in the 
Neumann case subsection. Following the same steps as in Section [5] and using (1551) 
we give the continuous spectrum contribution the form 



— / \(/ M DU,9u 



Sfc 



1 — cos(ap)q k e 1 ' 



R 2 



where 



X(k,p) 



1 - (q k e iak ) 2 
k 1+a p 1 -' T {2p + k) 



dkdp , 



(60) 



{p + k) 2 

We write this as the limit for e — > + of the integral restricted to k € (e, +oo) and 
expand the denominator into geometric power series (note that \q k \ < 1 for k > 0) 



1 - (q k e iak Y 



71=0 



'k.c 



iak\2n 



getting 



oo oo 



-L lim J ' M p [ X (k,p)Kz s k V qle ma 

ilf* £-5-0+ J J f— ' 



nG2N 



a fc cos(ap) J] g[V" ofc 



MG2N-1 



dkdp . 



We write the n = term separately as 

£ °o = ^3 / x(k,p)M p M k \s k \ 2 dkdp. 



(61) 



R 2 



For other terms we observe that 



N 

E 



qk& 



ik I 2?; 



< 



kl 2 M fc 2 

i-M 2 



< const Mt 



l + fc r 



m = 2 (D) 
to > 2 (N) 



(62) 



(see (|93|) and (J9TJ)) , so we can pull the infinite sum sign outside the integral by the 
Lebesgue dominated convergence theorem to obtain 



cont 

^ n 



^3 £ !™ + E / M *> / *( fc ' ^ (* e * 

ne2N l_e 

OO [ OC 

^ £ lim Y. J M pCOs(ap) J X (k, P )s k (q k e lak ) n dk 



dk + c.c. 



n£2N-l q 

We now split the analysis into separate cases 



dp 



dp . 



(63) 
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7.1 Dirichlet case (a = 1) 

In this case when (|5^|) is multiplied by x(k,p) the k~ 2 singularity in k = is 
canceled by k 2 from the function %. Therefore here the intermediate step with 
nonzero e is not needed and (p51) should be read with e = 0. Moreover, there is 
no bound state here, so this formula represents the total Casimir energy e a . If we 
write (na) 3 e mak — (~id k ) 3 e lnak and integrate three times by parts we find 

oo oc 

(na) 3 J X (k, p)s k (q k e lak ) n dk = -i^ - i J d 3 ( X (fc,p)«* tf) e mak dk , 
o o 

and then obtain (the first term on the r.h.s. above is imaginary and falls out) 



-a c oo 



67r 3 a 3 I 

I „£2N R ^ 



E ^3 / M v d l (*( fc > P> k «fc) ^"^ d P 



J2 ^ f M p cos(ap)d 3 k ( K x(k,p) S kq^)e mak dkdp\ + c.c. (64) 

■OlXl 1 ^ I 



n€2N-l R2 



Let now Q be the intersection of 3t+ with an arbitrary neighborhood of zero. We 
now use the results of Appendix [B] to infer that the following successive three 
operations on this formula lead only to the neglect of terms of order o(a~ 3 ): 

(i) replacement of dl (x(k,p)s k q k ) by s k q k dlx(k,p); 

(ii) restriction of the integration region to f2; 

(iii) replacement (in the restricted region) of M p s k q^ by MoSoSo = 1- 

In this way we arrive at 

2 ^ 1 fp 2 (k-3p) . . ,, „ , 
£ ° = £oo + ^ 2. -1 (fc+ )5 sm(nak)dkdp 

„£2N J Q V 1! 

-^3 E hj P {k\pf cos(ap)sin(nak)dkdp+ (1) . (65) 

nG2N-l ^ V - P/ 

The integrals are bounded uniformly with respect to a. It is now sufficient to show 
that they have well defined limits for a — > oo; then those limits determine the 
a~ 3 -term. 

Consider the integrals (with k,p > 0) 

f ( P 2 4p 3 A 
C(n,£,a) — / I 77 -7 — 7- 777 sin(nafc + £ap)dp dk , £ = 0, ±1 . 

7 \( fc +P) [k+P) J 
fc+p<l 
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It is easy to show that with the choice $7 = {k,p > 0,k+p < 1} the integral in the 
first line of (l65l is this integral with £ — 0, while the integral in the second line 
of (pj|) is one half of the sum of integrals with £ = ±1. By the change of variables 
r = a{k + p), t — p/(k + p) we bring this to the form 

1 a 

C(n, £,a)= I (t 2 - At 3 ) f - sin [r(n + t(£ - n))] dr dt , 
o o 

and find 



lim C(r 

a-^oo 


f-f for£ = 
l,e,a)= Wl 1 "' 3 " 4 "l for£- 


= 0,1, 

: -1. 


U 1/ ' (n+1) 3 (n+l) 4 J IOTi 


Using these results we get finally 


1 


[v 2 i v ( x 3n V 


+•(?) 




Lk» 8 Jk-A(»+i) 8 (»+DVj 




c(4) i ofM - " 2 i .r 1 


~"°° 167rV ' °U 3 J "« 1440a 3 1 ^ a 3 



where Eqq is defined in (1611) . 

7.2 Neumann case (a = —1) 

We now use formula (JS5")) with the replacement in notation q k e %ak = q k e lak , where 

■*0 ~ iafc 

Q- = — — - , a = a - a, g fc = q k e , 

ttMq 

which has the technical advantage that q' — 0. After this modification the general 
scheme is very similar to the Dirichlet case. Integration by parts gives expansion in 
1/a but at the end we shall translate it to the 1/a expansion. Integrating by parts 
we obtain boundary terms, for which in the present case the addition of c.c. terms 
must be taken into account before the limit e — > + is performed. For example, 
the first integration by parts in k in the first line in (1631) gives a term proportional 
to (before p-integration) 



lim + E lx^p)^[ Sk q n k {e)e m ^} 



«G2N _„ 

= - lim x ( e ,p) 3m [a e In (l - gf ~ 2 '<" ' - ' ' 



e 



where in the last equality we used the fact that in the neighborhood of zero 
I In |1 — g 2 e 2lae || < In cons '"' (see (HH)). After integrating by parts three times 
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in similar way we get the expression 



lim 

E^0 + 



oo 

£ / ' X {k lP ) Sk qle mhk dk + 
n£2wf 



np n 2 3JqC(3) 

'aM Ad, 2 M 2d 3 nM 2 



lim 



oo 
n62N V ; "I 



a fc ^ e mafc dfc 



c.c. 



For the sum over odd natural numbers the computations are similar. The integrals 
over p we treat in a similar way as in the Dirichlet case (steps (i) to (iii) , for their 
permissibility here see Appendix [Bj . The integrals of 9|x over the neighborhood 
of zero go similarly as in the Dirichlet case. In this way we find 



Cm 



- + 



6ir 2 M a 487rM a 2 



n-> 



1 /qC(3) 

6tt 2 8tt 4 M 2 



where 



C M = / P M p d P . 
o 
In the calculation we used the relations 



1440 



OO OO 

-Mprfp = q ' / P-^p cos(ap)dp = — r^ + °(^3 J 



<?)• 



(66) 



The first equality follows from normalization of /. We now return to the bound 
state contribution to the energy. From the second term of (|59|) . (14 3 p and (|44|) . 
using (j6"6")) we get 



abound 



Cm 



1 



67r 2 M a 6?r 2 a 3 



'(?) 



Adding the two contributions and changing the expansion parameter to a we 
obtain (|55J) . 

8 Local properties 

In this section we consider the local algebras of fields supported outside the regions 
of support of the potential V. For the initial (unsealed) models this means that 
the z-support of fields is outside the set (—b — R, —b + R)U(b — R,b + R) 7 but 
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for A — > eventually every support outside the planes z = ±b is admitted. Fields 
thus supported are also in the algebra of fields of the models with sharp boundary 
conditions at z = ±6, so one can also consider sharp boundary conditions for them. 

We recall from [T] that we use the initial value fields (smeared on a Cauchy 
surface of constant time) &(V), where V is a pair of real test functions (v, u), and 
X(u) = $(0, u), P(v) = $(v, 0) have the interpretation of canonical variables. For 
the present choice of the algebras the test functions are assumed to be in the space 
of smooth functions of compact support outside z = ±6, which we denote by XV 
The algebra of fields is formulated, more precisely, in Weyl form, which means 
that rather than $(V) elements W(F) = exp[i$(V)] are used, and the fields $ 
are defined on the level of specific representations. The states on the algebra are 
given as normalized positive linear functionals on the algebra of Weyl elements. 

With the Hamiltonians of the perpendicular motion h za and hf a as defined 
in (Unj) and flU} we denote 

h 2 a = h 2 za + hl, lh*] 2 = [hf a ] 2 + hl, (67) 

where — h\ is the two-dimensional Laplacian in the directions parallel to the plates. 
Then the ground states of the fields corresponding to the models proposed in the 
present article and to the sharp boundary conditions are given, respectively, by 

u a (W(V)) = exp [- \\\ 3a (V)\\ 2 } , ug{W(V)) = exp [- \\\fi [V)f] , (68) 

where 

j a {V) = hTv-ih-^u, j*(V) = [h») 1/2 v-i[hB]- 1/2 u. (69) 

(To be precise, to obtain formula (JTJ) for the Casimir energy we start in [5] with 
the directions parallel to the plates restricted to a box, whose size then tends to 
infinity. This may be shown to reproduce the states given above, but we omit this 
step here.) The ground states w a ,\ of the scaled models are defined analogously 
with the use of h a ,\- 

We show in this section that the scaled states reproduce in the weak limit 
the sharp boundary state: 

lhncj a , x {W(V))=^{W(V)). (70) 

A — HJ 

Also, we consider the limit of the local energy density. 

8.1 Local limit of states 

There is ||j a (^)|| 2 = (v, h a v) + (u, h^u), thus to prove ([70]) it is sufficient to show 
that ((p, h^\ip) tend to (ip, [h^] iji) for A — > and ip, ip G XV For such tp one has 
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for sufficiently small A. Therefore (<p, /i o ,aV0 = ~i { Pih a \/S.'iji) and similarly for 
ha , and the problem is reduced to the h~\-case. Moreover, ip and -0 are in the 
domain of h]_ , so we can write (ip, h~\tp) = (/i^ ' ip, h/ h~\h/ (hj_ ip)) 
and similarly for h^ . In this way the problem is reduced to the following: 



1/2, -1 ,1/2 , 1/2 r.fl-i -1,1/2 



lim^Xl^ =^ >?] *!' ■ ( 71 ) 



A-s-0 



To show this we first observe that \\h/ h \h/ | < 1, so it is sufficient to perform 

the limit between vectors from the total set of the form x{$) = X±(%±-)Xz{x 3 )- 
Using the spectral representation of h\ one has 

(X, hTKXhTp) = I (xz, , W±l P*) XTW)P1(P±) d 2 Pl _ . (72) 

The scalar product under the integral is bounded by ||Xz||||Pz|| an d for each p±, 
by the result of Section 02 tends to analogous expression with h z(l) \ replaced by 
hf a . This is sufficient to perform the limit, which ends the proof. 

8.2 Local energy density 

The state uj a (unlike the state w,f ) is defined on the whole algebra. In the language 
used in [T] the energy density in this state (properly normally ordered with respect 
to the vacuum) is given in the whole space by point-splitting procedure by 

£ a (x) = T a (x,x), (73) 

where T a (x,y) is the distribution defined by 

T a (^,V) = \(tp,{h a -h)i>) + \(W^,(h- 1 -h- 1 )^^), (74) 

with scalar product between the two gradients understood in the second term. 
The test functions are taken to be real. If one takes into account the translational 
symmetry in the directions parallel to the plates one realizes that the x±,y±- 
dependence of T a may be only through the difference x± — yj_. The removal of 
point-splitting in these directions means putting this variable equal to zero, or 
integrating the 2-dimensional Fourier transform of T a over all space of p±- variables 
of the spectral representation of h±. In this way one finds (now x, y are variables 
in the direction orthogonal to plates and (p, -0 are one-dimensional) 

£ a (x) = T za (x,x) , (75) 
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where 

+fl(v, m a +Pir 1/2 (ht+pir 1 ^) 

+ (V, [(hL +p1)~ 1/2 (hi + ft)- 1 ' 2 ]*') }dV , 

with prime in <p' ,ip' denoting the derivative. The integral is easily carried out 
explicitly and one finds 

r *«fa *) = ^(y, (?4 " ^)V0 - -^fe>', (ft,« - W) ■ (76) 

Our objective is to find the limit of the energy density in the scaled models. 
Thus following the introductory remarks of the present section we assume now the 
supports of ip, ip to be outside the set (—6 — R, —b + R) U (b — R,b + R) . For such 
functions there is h za ip = h z ip, so 

T za (^^) - -^-(y, (h ta ~ h,W) - ^-(<p', (h za - h z )4,') . (77) 

Z47T 57T 

Consider the general element (ip, (h za — h z )ij}) with test functions in the assumed 
class. For any non-negative real numbers a, b one has the identity 

oo 

a - b = -^J{wT^)-WT^)} r2dr - 



Using this and the spectral representations of h za and h z one finds 

(<p, (hza ~ h,W) = | / (^ [G(~r 2 ] - Go(-r 2 )]^j)r 2 dr . 
o 

To find the integral kernel [G(~r 2 ) - G Q (-r 2 )}(x,y) of [G(-r 2 ) - G (-r 2 )} one 
needs only to transform the kernel ([2~T]) to the position space: 

For x, y in the assumed region the Fourier integrals may be closed to contour 
integrals in the complex plane (the half-circular contributions vanish either in the 
upper or the lower half-plane due to the bound ()83[) ) and evaluated by residues. 
Because eventually we are interested in removal of point-splitting we assume that 
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x and y are in the same connected part of the region considered. We denote 
T(— r 2 ) = . I , where A and B are functions of r, and get 



[G{-r 2 ) G (-r 2 )](x,y) = -^-\f( ir )\ 2 x 



r i 

' e -\*+v\r[Acosh(ar) + B] for x,y > b + R 

or x,y < —b — R, 

e~ ar [A cosh ((« + y)r) + oB cosh ((x - y)r)] for x,y e (-6 + R, b - R) . 

We use this integral kernel in (|77|) . evaluate derivatives by parts and remove the 
point splitting. In this way we find in the assumed regions 

£ a (x) = T za (x,x) 

2e- 2 ^ r (Acosh(ar)+B)\ ? 3 for \x\ > b + R, 

* > \f(ir)\ r dr , 
e- ar (2Acosh(2xr)-aB)j tor \x\<b-R. 

We now want to consider the limit of this local energy in the scaled version of 
the model. Because of the appropriate convergence of the integral this limit may 
be performed inside the integral. The scaling of f(ir) and T(—r 2 ) follows from 
Section[3] As for the scaled model the excluded position set shrinks to ||x|— b\ < XR, 
in the limit all x ^ ±6 are admitted. A straightforward calculation yields 




-(2\x\—a)r 

l)r + 0e~ 2ar \ r?, dr j 

2ar ' 

With the help of Eq. (|94| in Appendix |A] we get 

a 



oo , g-(2|a:|-a) 

)im. £ at \{x) = ——2 / \ e (2x-a> + e -(2x+a) 
"> 1 _ g-2ar 



for |x| > b , 
for \x\ <b. 

for \x\ > b , 



167r 2 (|x|-6) 4 

£a 0) = lim £a,\{x) = { n 2 ^ a 

7 — > tt, 7T for \x\ < b . 

1440a 4 ^ 16TT 2 (nb-x) 4 ' ' 



A^0+ 



(78) 
where 2% + 1 denotes the set of odd whole numbers. 

Let us stress once more: the above result holds in the distributional sense 
only for test functions supported outside x — ±6, i.e. for functions in that class 
there is 

lim / £ a \(x)(p(x)dx — / £,f (x)ip(x)dx . 
>^o+ J : J 

For functions not in this class our assumptions leading to the above result cease 

to hold. Nevertheless, we shall attempt some comparison with our global results. 

For that purpose let us denote 

a 

167T 2 IE 4 



£ B {x) = - " , for x £ . (79) 
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Then ([751) may be written as 

£»{x) = £ B (x + b) + £ B (x -b)+ £ B (x) , (80) 



where 



for \x\ > b , 



167r 2 (|a;|+6) 4 

1440a 4 ^ 167r 2 (n6-x) 4 IOT Fl^°- 

n6(2Z+l)\{+l,-l} 

It is easy to see that £ B - mt {x) = a~ 4 F(x/a), where F(y) is an absolutely bounded 
and integrable function. Thus for large separation of plates the energy outside the 
plates is concentrated in the first two terms in ([50)1 . Therefore (|79p has the inter- 
pretation of the energy density produced around one single plate, while £ B int (x) 
may be regarded as the energy density (locally outside the boundaries) of the in- 
teraction. When integrated over x € K it gives 

7T 2 

^ int = "l440^ (82) 

for both Dirichlet and Neumann cases. 



9 Discussion 

The models considered in this article do not pretend to describe the details of the 
interaction of quantum field with macroscopic bodies in realistic way. Their merit 
comes from the fact that (i) they are consistent with the restrictions imposed by 
the general algebraic analysis of any quantum backreaction setting pQ while (ii) 
being simple enough to allow explicit calculations, and (iii) approximating sharp 
boundary conditions in a controllably localized (in physical space) way. The last 
point distinguishes them from the class of models constructed in [2] and allows the 
comparison of the results of the global and local analyzes. 

We summarize the results and lessons to be drawn from them. 

(i) The models discussed in the present article are defined on the Weyl algebra 
of the free theory, and the interaction introduced by the 'nonlocal boundaries' 
does not lead out of the vacuum representation of the theory. The energy of the 
field (as defined by the free theory) is finite in the ground state (as defined by 
the interaction with 'boundaries') without any arbitrary 'renormalization'. The 
change of this 'Casimir energy' with the variation of the position of the boundaries 
determines the backreaction force. 

(ii) The scaled interaction with the 'nonlocal boundaries' approaches in the scaling 
limit the sharp Dirichlet/Neumann conditions in the Hilbert-Schmidt norm sense 
for the resolvents. This implies norm convergence (or strong convergence) for any 
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continuous, vanishing in infinity (resp. any bounded) function of the first-quantized 
hamiltonian. The nonlocality of the boundaries is under control and tends to zero 
in the limit. 



(iii) The Casimir energy per area of the scaled models obeys the scaling law (|56|) . 
Thus it is governed by the a-dependence of the energy given in Eqs. ([57)) and (|58p. 
The Casimir force per area is minus the derivative of those formulas. We have not 
discussed this point, but one can show that one can differentiate terms in these 
formulas one by one with o(a~ 3 ) going over to o(a~ 4 ). Thus one finds 



480a 4 + 
1 


I s + 8tt 


Io 

2 M$a 4 


(3C(3) 

^ IT 2 


L ' 247rM A< 



det_ -480^ + °bJ' (D) 

~du~ 1 , 1 , J o (K{$) , ^ tt 2 /A 

480a 4 Va 4 



i)-7^-r + «(- ■ (N) 



One finds that in Dirichlet case the force has a well-defined limit, but in Neu- 
mann case depends on the model and diverges for sharp boundaries. This model- 
dependence occurs despite the fact that the models approximate in many respects 
the sharp boundaries very well (globally!). Neumann case models with Iq = 
(which are among those admitted by our assumptions - see Appendix C) have 
faster convergence property (see Eq. (|3"5j)') and for them the additional a~ 4 term 
in the force vanishes. 

(iv) Algebras of fields localized outside the interaction region (test functions with 
supports not intersecting with that region) admit free vacuum state, ground states 
of our models, as well as ground state of the sharp boundaries. Restricted to these 
algebras our ground states tend weakly to the sharp boundaries states in the 
scaling limit. 

(v) The local energy density is unambiguously defined in our models in the whole 
physical space by point splitting, with no ad hoc later renormalization. When 
restricted to regions not intersecting with the interaction area (which in scaling 
limit just means not intersecting with the boundaries) the local Casimir energy 
has a well defined limit given by a smooth function Eq. (l78l) . This limit density 
has universal, mo del- independent form. Thus the 'bulk' contribution to the total 
Casimir energy has this universality and the mo del- dependent terms in Neumann 
case turn out to be squeezed in the limit inside the boundaries. The hope that 
these non-regular contributions may be removed in a model-independent way is 
therefore not justified. Further confirmation of our interpretation supply formulas 
and remarks ending the last section. 

Local energy density has been discussed by various authors before @]-[B], but 
usually with the use of some regularizations, often of not quite clear status. The 
results to be found in literature are not consistent. The Dirichlet case for the region 
between the plates is discussed in [5] and modulo some infinite renormalization 
agrees with ours (formula (2.32) in that reference). On the other hand the authors 
of [6] obtain a different result (by a rather indirect way of 'regularization' and 
removal of cut-off). We are not aware of a complete rigorous discussion resulting 
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in our formulas (175)) - (152)) . It is also worth noting that the density £,f (x) can be 
also obtained directly by the use of sharp boundary conditions hamiltonian hf a 
instead of h za in (fTTj) . This amounts to the use of the difference of (|55)) and (|52)) in 
the calculation. However, let us stress once more, this limit value of the density is 
correct only if smeared with a test function with support not touching the borders. 

We are grateful to the referee for careful reading and constructive editorial remarks. 

Appendices 

If not stated otherwise we work here with the same assumptions as stated in (|lip 
and (HI. 



A Integrals and estimates 

We gather here a few technical, separate points used in the main text. 

(i) We recall that / is a smooth function supported in (— R, R) if, and only if, 
its Fourier transform / is an entire analytic function satisfying the estimates 

~ ,. const(iV)e fl l : ' mu l 

for all u e € and N 6 N. 

Therefore M p , which is the product of the Fourier transforms of f(z) and 
of f(—z) = f(z), extends to the analytic function M u on (D satisfying similar 
estimates with R replaced by 2R. The inverse transform M of M p is smooth and 
supported in (— 2R, 2R). Using this we find for a > 2R 

cos(ap)M p dp = J^ [M(a) + M(-a)] = . (84) 

As M u is even, — | — — also extends to an analytic function and with the use 

of estimates on M u one finds by closing the contour of integration (as always for 
a > 2R) that 

cosMCM-M.) 

pZ _ fcZ 



2r, 



(ii) In order to prove the estimate (|50)l . we start with Dirichlet case, we note 
that (remember that /'(0) = 0, as / is even) 



f(Xp) - /(o) 



r(pe)(A-e)« 



<Ay |/"(pe)|de. 



Moreover we have (norms of functions of p as elements of L 2 ) 



\f"(pSM 



< / ||/"(K)||||/"(K)||^^ 

<0,A) 2 



< 11/'' 



(C|)" 1/2 rf^<constA, 



(0,A> 2 



P 



2 2 



where in the first step we have used the Schwartz inequality. Now because 

is bounded we get the estimate for Dirichlet. The same considerations show (lUll 
The proof of Neumann case in (J30I) is almost the same but with the use of the 
formula 

A 

f(Xp)-f(0)=p (?(pZ)di;- 



(iii) For the function Ik, (jMJ), we first observe that for some < eo < 1 there is 

Jo < 2nRM (l - e ) . (86) 

This follows from the identity 

oo 

sin 2 (Rp)M p 



J 



-dp = 2ttRM - I , 



(used already in the normalization of the bound state, see (J43j) , (|44j) ) , as the l.h.s. 
is strictly greater than zero. Further, we need the following estimates: 

const (fc*) 

(87) 



h> 



h< 



k 2 
const 



k> 0, 



k "= (1 + k) 2 ' 

with arbitrary k* > 0. We write Ik as a principal value distribution calculated on 
test function M. In position space, using evenness of /, we have 

oo 

I k = V2^ [xM(xf-^±dx. 
J kx 



2(i 



Integrating once by parts we get 



M(0)+ / M'(x)cos(kx)dx 
o 



(89) 



The estimation from above is now trivial, whereas for the estimation from below 
we use the Riemann-Lebesgue lemma, the assumption that I k > for k ^ (see 
<j37j)) and continuity of Ik- Expanding further in powers of 1/fc the integral in (|89| 
we find moreover 

(iv) The Sfc function, defined in the end of Section @J is smooth for k > and 
satisfies 

|9fcSfe| < const(n) (1 + k)- {n+ ^ . (91) 

To show this we note first that \M k — iNk\~ l is bounded in a neighborhood of 
k = as Mo > 0. Outside this neighborhood we have \M k - iN^ 1 < \Nk\~ 1 < 
const (1 + k)~ a due to (|38|) and the first bound in (|87|) . On the other hand due 
to d90]) there is \d%N k \ < const (1 + £;)~ ( ' l+1) for n > 1 except for n = 1 in the 
Dirichlet case, when this is replaced by \dkN k \ < const. As M& is Schwartz, this 
ends the proof. 

It now follows that q k — s^M^ and all its derivatives are bounded by Schwartz 
functions. 
(v) For both (D and N) cases and k > we note the bound 

1 , N 1 + k , s 
— < const (a) — — . (92) 



\l-(q k e iak ) 

At k — there is 1 — (go) 2 = but the first derivative of the denominator at 
zero does not vanish (in Neumann case use (j5r?)) ). which is sufficient for (|92"|) in a 
neighborhood of zero. Outside that neighborhood, mainly due to (|57|) . it follows 
that 

1 1 M fe + ^fe ^ 

1 ^7 < i i fo = \r9 < C0I1St I 

|l-(%e» fe ) 2 | " 1-1* I 2 A^l 

which ends the proof. We also note that for k > we have 

< const i__ \'"~~' Z- ( !): '») 



i - 1* 




where the Neumann case depends on the behavior of Ik at zero (m = 2 for Io =/= 
and m = 2 + 2r when Ik — k 2r , r > 1 as /& is even). 
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(vi) Finally, we note the following identity using a known integral representation 
of the Hurwitz zcta function: for a > there is 



1 i re_^_ d J_ (■( 4 « \ = Y (a + 2an)- i , (94) 

o n -° 

which is needed for the calculation of (j78j) . 

B Operations (i) — (iii) from Section [7] 

In this appendix we prove the admissibility of the three operations (i) - (iii) per- 
formed in Section 17.11 for the Dirichlet case and mentioned in Section 17.21 for the 
Neumann case. The key tool for this is the following simple lemma: 

Let c„ for n G JV C N be complex measurable functions on D C R. 
If Ylme.A' l c ™(^)l ls hitegrable on D then 

lim V / c n (k)e inak dk = 0. (95) 

a— *-oo -^ — ' / 

The proof is straightforward and uses the Lebesgue dominated convergence the- 
orem and the Riemann-Lebesgue lemma. Throughout this appendix S^p. denotes 
some Schwartz function in fc variable (each time it may be a different function). 
All we have to do is to check if the assumption of the mentioned lemma is fulfilled 
for the appropriate expressions. 

B.l Dirichlet case 

Let us first consider the following part of the expression (|64|) 

71G2N 



5^ -T / M P d l (x(*, PH <&) e mak dk dp ; (96) 

^ /- O TXT ** 



HiJ. 



we shall need its real part multiplied by — 1/67T 3 in (|64[) . 

(i) To justify the validity of operation (i) of Section 17.11 it is sufficient to check 

whether the sequence of functions 



c n {k) = n- 3 / M p dl{s k q%x{k,p)) - Sk q^d 3 kX (k,p) 



dp 



satisfies the assumption of lemma (|95|) . For this we note that using properties of 
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n" 3 



Sk and qk described in Appendix A (iv) we have 
dl{s k ql X {k,p)) ~ s kq tdl X (k lP ) 

2 

< ^(n^lgfcl"- 1 +n- 2 \q k \ n - 1 ) ^ |^x(*,P)l 

1 

^ fe n- 2 (n - l)|%| n - 2 Y, \ d lx(k,p)\ + y k n- 2 (n - l)(n - 2)| % r 3 | X (fc,p)| . 

3=0 



(97) 



A simple calculation yields 



\dlx(k,p)\< const £-i, i = 0,l,2, (98) 

fc +p 

and therefore 

M p |d£x(fc,p)|dp < const A: 2 - J log(l + A;- 1 ) , j = 0,l,2. (99) 

o 

Using this and remembering that \qk\ < 1 we obtain 

\c n (k)\ < y k log (1 + fc- 1 ) (n" 3 + n" 2 + ^fc^r 2 + fc 2 |g fc r 3 ) , (100) 

where the third and the fourth terms inside the parentheses appear for n > 2 and 
n > 3 respectively. The summation gives 



J2 |c„(fc)|<^log(l + fc- 1 ) 



n€2N 



1 + fc log ' 



i-kfcIV i-kfel 5 



(101) 



The estimate (p?3")) shows that the integral of the r.h.s. over R + is finite, thus the 

assumption of lemma (1^51) is satisfied. 

(ii) For the operation (ii) the assumption of the lemma is fulfilled because 

\M pSk q^dlx(k, P )\ < &P [^t M p^ e ^K\«) > 

with f2 as defined in Section I7TT1 

(iii) The operation (iii) is admissible because 

\MpSkqk ~ ^o^oQo I < const p + const (1 + n)k + const (1 + n)pk 

and kdlx(k,p),pdlx(k,p) are both integrable on f2, so the lemma holds also in 
this case. This ends the proof for the terms (f96|) . 
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For the expression with the sum over odd natural numbers and with cos(ap) 
we use the following modification of the lemma. Let c n for n <E JV C N be complex 
measurable functions on D C E 2 . If ^2 n€ ^y |c n (fc,p)| is integrable on D then 

lim V [ c n (k,p)e mak e ±lap dkdp = 0. (102) 

a— »oo -^ — ' / 

Now, almost the same considerations as before show the admissibility of operations 
(i) - (iii) for this part of energy. 

B.2 Neumann case 

First we note that for the terms which we consider here, the limit over e can 
be easily performed. This follows from the estimations below. Therefore we use 
the same lemma as for Dirichlet case, i.e. (|95|) and (|102[) . but we recall that we 
now replace q^e lak = q^e lak , as mentioned at the beginning of Section ITT2"1 With 
this modification the estimation (|5T|) for the terms in the sum (|M|) is still valid, 
but the bounds (|55|) change: k 2 ~i is replaced by p 2 ~ J . In consequence there is no 
k 2 ~ j factor in front of log in ([55]) and no k, k 2 factors in ([TOD]) and (fTUT|) . With 
this modification the sum of n~ 3 , rT 2 and rT 1 terms in (|100|) is still sufficiently 
well bounded, but the sum of n° terms is to singular (no k 2 in the last term in 
(jl01() ) . Therefore these terms need a more detailed treatment. In fact, they can be 
estimated by const |sfc||gfc| ra ~ 3 |<?£.| 3 log(l + A: -1 ) (prime denotes here the derivative 
with respect to k) and their sum over n is bounded by 

const l^f wi + fc- 1 ). 

This function is indeed in L 1 (1R + ), since outside the neighborhood of zero, using 

\sk\\q' k \ 3 

1-14 1 
is enough for I ^ case (see (|93|) and the comment after it) and if Ik — k 2r , 

k— >0 

W I 3 

r > 1 (ife is even), then — -tt ~ k 2r ~ 2 . The discussion of admissibility of 

1 - \qk\ fe ^° 
operations (ii) and (iii) goes in the same way as for Dirichlet case (estimates hold 

also for Neumann case, with qk instead of qk). The analysis of the expression 

with the sum over odd natural numbers and with cos(ap) is also analogous to the 

Dirichlet case. 



C Comments on the assumed class of models 

The models discussed in this paper are based on functions / subject to the condi- 
tions formulated in (IPTl) . (|T2"j) . (|57|) and (f55|) . In this appendix we exhibit a class of 
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93|) . we have k „ < S*k, whereas in the neighborhood we have q' = which 



functions conforming to them. It is sufficient to satisfy (jTTJ) and (J3TJ) as the other 
two are then achieved by simple rescaling of the function by a constant factor. 

First, we note that each even function with the assumed support which in 
addition is real, non-negative and monotonically (weakly) decreasing for positive 
arguments satisfies the demands. Indeed, for each function the last condition in 
(fTTj) is fulfilled. Moreover, it is easy to see that then M{x) (being the convolution of 
the function with itself) is also even, positive, compactly supported and decreasing 
for x > 0. Thus from (|55|) . since for k ^ there is 



< - / M'(x)dx = M (0) 



M'(x) cos(kx)dx 



we have Ik > for k ^ 0, which ends the proof of (13"T|) . 

For each function in the class defined in the previous paragraph there is 
Jo > 0, which is due to the positivity of M(x) (see (|55])). We now extend the class 
to include also functions for which Iq = (see Discussion). Let / be a function in 
the class of the last paragraph and define a new function 

f{z) = f{z)-^f(z-r)+f{z + r)) 1 
where [i > and r > R > 0. One finds 

M r p = \f r (p)\ 2 = (1 - 2/icosM) 2 M p . 

Using this and taking into account (|55|) one has 

f Ml -MI / „\ sin(rfc) 

J k = I 2 _ k 2 d P ={ l + 2 M J4 - 4^/i— |p (1 - M cos(rfc)) M k . 

We impose the condition Iq = 0, which is a quadratic equation for /z: 
2(7 + 27rrA/oV 2 - 47rrAf 0) u + I = . 

For sufficiently large r the equation has two roots, and we take the smaller one, 
which is less then 1/2 and for large r tends to zero. Then 

II = 47rr/i(l - M )Mo^ [1 - »7(rfc)£(&)] . 

where 

sin(rfc) l-/icos(rfc) I Af fe 

77(rfc) = , 4(fc) = . 

rk 1 — fi MqIu 

It is an exercise in function analysis to show that for sufficiently small /i there 
is r)(u) < 1 for all a/0. Below we show that / may be chosen such that also 
£(fc) < 1 for k ^ 0, and then I r k > for k > 0, which is the condition ([57)1, To 
reduce the size of the support of f r one can use the scaling defined in Section [31 
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Consider now the non-negative and even function £(fc) for positive arguments. 
Suppose that £"(0) < (see below). Then there is ko > such that for k € (0, ko) 
it is £(fc) < 1. For fc > fco, using ([571) . we have £(fc) < const = £ max . We now need 
to improve the estimation of r\ for k > fco- It is easy to see that if r has been chosen 
large enough then r](rk) < -r^— for k > k . 

Finally, we have to choose / so as to satisfy £"(0) < 0. Let / be defined as 
f(z) = fa = const on (-R, R) and zero outside. Straightforward calculation gives 

The function used here needs 'rounding the corners' to be in the class of the second 
paragraph. But this may be made by a small local variation, and both M& and Ik 
and their derivatives depend continuously on such small variations of / (for Ik see 
which is sufficient to conclude the proof. 
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